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ABSTRACT 

Context. The Lupus I cloud is found between the Upper-Scorpius (USco) and the Upper-Centaurus-Lupus (UCL) sub-groups of the 
Scorpius-Centaurus OB-association, where the expanding USco HI shell appears to interact with a bubble currently driven by the 
winds of the remaining B-stars of UCL. 

Aims. We want to study how collisions of large-scale interstellar gas flows form and influence new dense clouds in the ISM. 
Methods. We performed LAB OCA continuum sub-mm observations of Lupus I that provide for the first time a direct view of the 
densest, coldest cloud clumps and cores at high angular resolution. We complemented those by Herschel and Planck data from which 
we constructed column density and temperature maps. From the Herschel and LAB OCA column density maps we calculated PDFs 
to characterize the density structure of the cloud. 

Results. The northern part of Lupus I is found to have on average lower densities and higher temperatures as well as no active 
star formation. The center-south part harbors dozens of pre-stellar cores where density and temperature reach their maximum and 
minimum, respectively. Our analysis of the column density PDFs from the Herschel data show double peak profiles for all parts of the 
cloud which we attribute to an external compression. In those parts with active star formation, the PDF shows a power-law tail at high 
densities. The PDFs we calculated from our LABOCA data trace the denser parts of the cloud showing one peak and a power-law tail. 
With LABOCA we And 15 cores with masses between 0.07 and 1.71 M© and a total mass of ^ 8 M©. The total gas and dust mass of 
the cloud is ^ 164 M© and hence ~ 5% of the mass is in cores. From the Herschel and Planck data we And a total mass of 174 M© 
and ^171 M©, respectively. 

Conclusions. The position, orientation and elongated shape of Lupus I, the double peak PDFs and the population of pre-stellar and 
protostellar cores could be explained by the large-scale compression from the advancing USco HI shell and the UCL wind bubble. 

Key words. Stars: formation - Stars: protostars - ISM: bubbles - ISM: clouds - ISM: dust, extinction - individual objects: Lupus I 


1. Introduction 


In the current picture of the dynamic interstellar medium (ISM), 
molecular cloud formation is attribute d to collisions of l arge- 
scale flows in the ISM (see review by Dobbs et"ar||2014[ and 
references within). Such flows can be driven by stellar feed¬ 
back processes (e.g. UV-radiation and winds) and supemovae. 
At the interface of the colliding flows, compression, cooling, 
and fragmentation of the diffuse atomic medium produces cold 
sheets and filaments that later may become molecular and self- 
gravitating and domin ate the appearance of the ISM as observed 
today (see review by |Andre et al. 2014). In this picture the 
fast formation (and dispersion) of molecular clouds and the of- 


^ The Atacama Pathfinder Experiment (APEX) is a collabora¬ 
tion between the Max-Planck-Institut fr Radioastronomie (MPIfR), 
the European Southern Observatory (ESO), and the Onsala Space 
Observatory (OSO). 

** Herschel is an ESA space observatory with science instruments 
provided by European-led Principal Investigator consortia and with im¬ 
portant participation from NASA. 


ten simultaneous onset of star formation within (see [Hartmann | 


et al.|2QQT||Vazquez-Semadeni et al.|2007[|Banerjee et al.|2009[ 

Gomez & Vazquez-Semadeni|2014| ) appears plausible. 

One example of such a large-scale flow is an expanding shell 
or super-shell around e.g. an OB-association or in general driven 
by multipl e stellar feedbac k of a star cluster or association (see 
review by |Daws^ 2013). Molecular clouds may then either 
form inside the wall of such a shell ( [Dawson et al.||2011| ) or at 
the interface region when two such shells collide with each other. 
The latter has been recently investigated by jPawson et al. (2015) 
for a young giant molecular cloud (GMC) at the interface of two 
colliding super-shells. Erom the comparison of CO observations 
with high-resolution 3D hydrodynamical simulations they found 
that the GMC assembled into its current form by the action of 
the shells. 

The Scorpius-Centaurus QB-association (|Blaauw|p^64} jde 


Zeeuw et al.|1999[|Preibisch et al.|2002{|Preibisch & Mamajek 

2008[ Sco-Cen) is the closest site of recent massive star forma¬ 
tion to us and it consists of three sub-groups with di fferent ages 
and well known stellar populations down to 2Mq (de Bruijne 
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1999| ). The oldest one is the the Upper Centaurus-Lupus (UCL) 
sub-group with an age of ~ 17 Myr harboring 66 B-stars. With 
an age of ~ 15 Myr the Lower Centaurus Crux (LCC) sub-group 
is somewhat younger and contains 42 B-stars. The youngest sub¬ 
group is Upper-Scorpius (USco) with an age of ~ 5 Myr and 
consisting of 49 B-stars. 

The feedback of the numerous massive stars in Sco-Cen 
probably cleared the inner region of the association from diffuse 
matter creating expanding loop-l ike HI shells around each of 
the sub-groups of the association ( [de Geus|1992| ). At the edge of 
the USco shell several dense molecular clouds with very young 
(< 1 - 2 Myr) stellar populations are found. Of those the most 
prominent ones are the Lupus I cloud (near the western edge of 
the shell) and the p Oph cloud (near the eastern edge). 

The Lupus I molecular cloud complex (for an overview see 
Comer6n|2008| ) is found at a distance of 150 pc and consists of 
a 2.6° X 0.6° (corresponds to a physical size of 6.8 x 1.6 pc) 
main filament extending in a north to south direction (Galactic 
coordinates) and a ring-like structure of ~ 0.6° in diameter wes t 
of the m ain filament (towards U CL; see [Tothill et al.||2009| ). 
Recently [Matthews et al. ( 2014| ) noted also two smaller sec¬ 
ondary filaments of which one is about half a degree long and 
runs perpendicular to the main filament and seems to connect 
with it in the south. The other one is about a degree long and lies 
south-west of the main filament extending from the southern end 
of the main filament to the ring-like structure. 

In this work we concentrate our analysis only on the main 
filament commonly seen in all observations. We will refer to it 
as Lupus I or the Lupus I filament. 

The cloud is found between the USco and the UCL sub¬ 
groups at a location where the expanding USco shell appears 
to interact with a bubble currently powered by the winds of the 
remaining B-stars of UCL. With its close distance Lupus I repre¬ 
sents a good candidate where we can study how such a collision 
process may form and influence new dense clouds in the ISM. 

Lupus I has been mapped as part of several large surveys 
like the Herschel Gould-belt survey ( Andre et al.||2Q10} jRyglj 
et al.|[MT3] ) and the Spitzer Legacy Pr ogram ’From molecu¬ 


lar clouds to plane t-forming disks’ (c2d; Chapman et al.||2007[ 
Merin et aL]|2008| ). These near-infrared to far-infrared surveys 


revealed the population of young stellar objects (YSOs) within 
the cloud showing that it is dominated by pre-stellar and proto- 
stellar cores indicating an on-going star formation event. jRyglj 
jet al.j ( |MT^ found that the star fo rmation rate (SFR) is increas¬ 
ing over the past 0.5 - 1.5 Myr and |Merm et al.j ( [2008) ) estimated 
a SFR of 4.3 M© Myr“^ for Lupus I from their Spitzer data. 

Extinction maps of Lupus I have been created by various au¬ 
thors using different methods. Cambres3^ ( jl999j ) created an ex¬ 
tinction map based on optical star counts with a resolution of 
a few arc-minutes. Using 2MASS data [Chapman et al.| ( j2007j ) 
created a visual extinction map with a 2' resolution. Also from 
2MASS data a wide field extinction map of Lupus I has been pre¬ 
sented by jLombardi et al.j ( 2Q08 ). It allowed extinction m easure¬ 
ments down to Ak = 0.05 mag, but had a resolution of S'. jMermj 
et al.j ( j2008j ) created extinction maps from their Spitzer data by 


estimating the visual extinction towards each source classified 
as a background star, based on their spectral energy distribution 
(SED) from 1.25 to 24 pm. Their maps had a resolution of 2 - 5'. 
[Rygl et alT ( j2013j ) created a column density map of dust emission 
in Lupus I from their Herschel data by making a modified black 
body fit for each pixel in the four bands from 160 - 500 pm. The 
resulting resolution was 36". 

Here we present a far-infrared and sub-mm analysis of the 
Lupus I cloud based on newly obtained APEX/LAB OCA sub- 


mm continuum data at 870 yum and complementing those with 
Herschel and Planck archival data. In Sec. we briefly describe 
the observations used and their data reduction process. Sec. 
gives an overview of the methods used to create column density 
maps and probability distribution functions (PDFs) from the ob¬ 
servations. We present and discuss our results in Sec.|^ In Sec.|^ 
we consider the surroundings of Lupus I and the influence of 
large-scale processes on the cloud. Finally, we summarize and 
conclude our paper in Sec.|^ 

2. Observations and data reduction 

2.1. LABOCAdata 

The sub-mm continuum observations of Lupus I were per¬ 
formed with the APEX 12-m teles cope located in the Chilean 
Atacama desert ( jCiisten et al.|2006j). We used the LA rge APEX 
Bolometer CAmera (LABQCA, jSiringo et al.|2009) ) which op¬ 
erates in the atmospheric window at 870yum (345 GHz). The an¬ 
gular resolution is 19.2"(HPBW), and the total field of view is 
11.4'. At the distance of Lupus I (150 pc) the angular resolution 
corresponds to a spatial scale of ~ 2800 AU (~ 0.01 pc). This is 
sufficient to resolve the structure of molecular cores. 

The LAB OCA observations of Lupus I were obtained in 
Max-Planck and ESO Periods 91 and 92, on 24th March, 31st 
August, 12th, 13th September 2013, and 31st March 2014 (PI: 
B. Gaczkowski). They were performed in the on-the-fiy (OTF) 
mode, scanning perpendicular as well as along the filament’s ma¬ 
jor axis with a random position angle to the axis for each scan to 
reduce striping effects and improve the sampling. The total ob¬ 
serving time was 11.3 hours and the weather conditions good to 
average. 

Data reduction included standard steps for sub-millimete r 
bolometer data, using the BoA software package ( jSchullerj2012j ). 
First, data were converted from instrumental count values to 
a Jansky scale using a standard conversion factor, then a fiat- 
field correction (derived from scans of bright, compact sources) 
was applied. Corrections for atmospheric opacity were derived 
from sky dips taken every 1-2 hours, and finally residual correc¬ 
tion factors were determined from observations of planets and 
secondary calibrator sources. Data at the turning points of the 
map were flagged, as well as spikes. The data from the indi¬ 
vidual bolometers were then corrected for slow amplitude drifts 
due to instrumental effects and atmosphere by subtracting low- 
order polynomials from the time-stream data. Short timescale 
sky brightness variations (sky-noise) was corrected by removing 
the correlated (over a large number of bolometers) signal in an 
iterative fashion. Then the time-stream data were converted into 
sky-brightness maps for each scan, and finally all scans com¬ 
bined into one map. 

The removal of the correlated sky variations is known to filter 
out also astronomical emission from extended sources, usually 
leading to negative artifacts surrounding bright emission struc¬ 
tures. In order to recover some of the extended emission, an it¬ 
erative source modeling procedure was applied, using the result 
from the previous iteration to construct an input model for the 
following iteration. To construct the model, the previous map 
was smoothed, all pixels below a pixel value of 0 set to 0, and 
the map smoothed again. In the following iteration, the model is 
subtracted from the time-stream data before de-spiking, baseline 
subtraction, and sky-noise removal, and added back to the data 
stream afterwards, before the new map is created. This procedure 
effectively injects artificial flux into the mapped area and is prone 
to create runaway high surface brightness areas especially in re- 
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gions without significant astronomical signal and low signal-to- 
noise. For that reason, the model image was set entirely to 0 in 
areas of the map with poor coverage, below a certain threshold 
in rms (basically the map edges), and the procedure was stopped 
after 25 iterations. Typically, structures on scales of the order of 
3-4' can be recovered, larg er scales get more and more filtered 
out (see simulations done by |Belloche et al.|2011| ). The gridding 
was done with a cell size of 6.1" and the map was smoothed 
with a Gaussian kernel of size 9" (FWHM). The angular resolu¬ 
tion of the final map is 21.2" (HPBW) and the rms noise level is 
23 mJy/21.2"-beam (hereafter the notation per beam in context 
of LABOCA will mean per 21.2"-beam). The resulting map is 
shown in Fig.[2 

In order to determine the total sub-mm flux of the entire 
LABOCA map, we integrated all pixel values above the 3cr noise 
level of 69 mJy/beam. This yielded a total flux value of 476 Jy. 
Our peak intensity in the map is = 1.37 Jy/beam. 



Jy/21.2"- 
beam 

0.42 

0.31 

0.21 

0.13 

0.06 

0.002 

-0.043 

-0.075 

-0.094 


Fig. 1. Lupus I LABOCA 870yum map with square-root inten¬ 
sity scaling in units of Jy/21.2"-beam. The positions of the 15 
cores that were found by Clumpfind are marked with white 
crosses. The white polygon marks the common area within 
which the mass of the cloud was calculated fro m the LABOCA 
and Herschel column density maps (see Section [4^. 


2.2. Herschel archival data 


The Lupus I cloud c omplex was observ ed by the Herschel far- 
infrared observatory ( [Pilbratt et al. 2010|) in January 2011 as par t 
of the Gould-belt survey ( [Andre et al.||2010[ |Rygl et al.||2013 ). 
The Photodetector A rray Camera and Spectrometer (PACS, 
Poglitsch et al.|2010 ) and the Spectra l and Photometric Imaging 
Receiver (SPIRE, Griffin et al.|2010| ) were used to map an area 
of 2° X 2.3°. We retrieved these data from the Herschel sci¬ 
ence archive and reduced them using HIPE vl2.1 ( |Ott 2010) 
for the calibration and the deglitching of the Level 0 PACS and 
SPIRE data. The Level 1 data of both instruments were then used 
to produce the final maps with the Scanamorphos package v24 
( |Roussel|2013| /parallel and /galactic options were switched on 
in both cases to preserve extended structures in complex, bright 
Galactic fields). The pixel-sizes for the five maps at 70,160, 250, 
350 and 500yum were chosen as 3.2", 4.5", 6", 8", and 11.5", 
respectively. 

The absolute calibration of the Herschel data was d one fol¬ 
lowing the approach described in [Bernard et al. ( 2010| ). Using 
the Planck and the IRAS data of the same field as the Herschel 
data one calculates the expected fiuxes to be observed by 
Herschel at each band and from that computes the zero-level 
offsets (see Table [^. Those are then added to the Herschel maps 
to create the final mosaics. 


Table 1. Zero-level offsets for the Herschel maps obtained from 
the cross-correlation with Planck and IRAS. 


Band 

[y^m] 

Offset 
MJy sr"^ 

70 

5.68954 

160 

14.3011 

250 

8.11490 

350 

4.06405 

500 

1.68754 


Eor the dust properties and temperatures considered in our 
analysis the color corrections for both PACS and SPIRE are of 
order 1-2% and hence negligible. 

2.3. Planck archival data 

We retrieved the all-sky maps of the High-Erequency-Receiver 
(HEI) in the 353, 545, and 857 GHz bands (corresponding to the 
EIR and sub-mm wavelengths of 850, 550, and 350/im) from 
the Planck legacy archive (release PRl 21.03.2013). Erom these 
data cubes in HEALPi^D format we extracted for each wave¬ 
length separately the area covering the Lupus I region and cre¬ 
ated gnomonic projected maps. The resolution in each of the 
three bands is 5' and the pixel size was chosen as 1.7'. 


3. Data analysis 

3.1. Column density and temperature maps of Lupus I 

Prom our three data sets of LABOCA, Herschel, and Planck we 
calculated column density maps to characterize the cloud in a 
multi-wavelength approach. Additionally, we constructed a tem¬ 
perature map from the Herschel data. 


^ http://healpix.jpi.nasa.gov/html/intro.htm 
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3.1.1. Column density and temperature map from Herschel 
SED fit with all SPIRE bands 


The standard and often practiced way to derive column density 
and temperature maps from the Herschel data is to fit an SED 
to the observed fl uxes of the Herschel bands for each pixel of 
the m aps (see e.g.|Alves de Oliveira et al.|2Q14[ [Lombardi et al. 


2Q14t Battersby et al.||2Q14| ). Here we fit the SEP to the three 

SPIRE bands 250, 350, and 500yum. We do not include the two 
PACS 70 and 160yum bands, because those observations were 
corrupted by stray moonlight and hence are not reliable for an 
analysis of large-scale structures. 

Because Lupus I is optically thin to the dust emission at the 
considered densities and wavelengths we can model its emission 
as a modified black body. We assume that the long wavelength 
emission (T > 250yum) of a pixel [i,j] comes from a unique 
species of grains being all at the same equilibrium temperature, 
and having a power-law wavelength dependent opacity. If Ly\X) 
is the monochromatic luminosity of pixel [i, j] at wavelength T, 
then it can be expressed as 


\X) = X 4^) (X AnBM, 


( 1 ) 


thin ice mantles and p = 10^ cm“^). Hence, our chosen value of 
/c(/lo) might still vary by a factor of about 2. 

The statistical error on both the final column density map 
and the temperature map consists of errors of the calibration, 
the photometry, and the SED fitting process. We conservatively 
estimated the sum of these uncertainties to be ~ 20 % for both 
maps. 


3.1.2. Column density map from Herschel SPIRE 250pm 
map 

To obtain a column density map with the resolution of the SPIRE 
250 pm ban d (i.e. EWHM 250 = 18''), we followed the technique 
described in Juvela et ST] ( |2013| ). Using the intensity map of the 
SPIRE 250 pm band and the previously calculated temperature 
map (at the lower resolution) to compute the Planck function 
it is possible to gain another factor of two in resolution 
compared to the SED fitting (to all three SPIRE bands) case. In 
this way the column density simply is 


N\ 


■H 250 


250yum 


H 2 


5y(250yum, r'’-') /C250/imyUH2 mn 


(3) 


where k is the dust opacity, p the emissivity index, By{A, T^) the 
black body spectral flux density for a dust temperature T^. The 
two free parameters and are the dust mass and tem¬ 
perature per pixel of the material along the line-of-sight. Here 
we adopt a typical a:(/Io) = 5.91 cm^g“^ with To = 350pm for 
dust grains with thin ice mantels and gas densities < 10 ^ cm“^ 
( jOssenkopf & Henning||1994| ). Using the available information 
about p from the Planck data, we found that the emissivity in¬ 
dex within the Lupus I cloud lies between ~ 1.6 and ~ 1.7. 
Therefore we fixed its value to yS = 1.65 as an average value 
within the cloud. Considering the low resolution of the p map of 
30' this is a reasonable approximation. In this way we also limit 
the number of free parameters in the fit making it more stable. 

After convolving the 250 and 350pm maps to the resolution 
of the 500pm band using the kernels from jAniano et al.| ( |2011 ), 
the modified black body fit was performed pixel-by-pixel. Erom 
the dust mass in each pixel the total column density for both 
dust and gas was then calculated as 


yyHsED ^ 
H 2 


MujR 

PU 2 


( 2 ) 


where R = 100 is the gas-to-dust mass ratio, ph 2 = 2.8 the 
molecular weight per hydrogen molecule, and mu the hydrogen 
atom mass. The resulting dust temperature values from the 
fit a t each pixel give the temperature map of Lupus I (shown in 
Eig.| 2 (a)|) at the resolution of the 500pm band, i.e. EWHM500 = 
36". 


Since the composition of the dust grains and their density 
is unknown, the choice of a particular dust model and thus a 
specific opacity law introduces an uncertainty of the dust opac¬ 
ity and thus of the resulting column density. In order to esti¬ 
mate it, we took different k(Ao) values at 350pm for grains with 
and without ice mantles and the three initial g as densities of 


10 ^, 10 ^, and 10 ^ cm ^ in the online table of 


Ossenkopf& 


Henning| ( |1994| ). They vary between 3.64 cm^ g ^ (MRN^with- 
out ice mantles and p < 10^ cm“^) and 11.3 cm^ g“^ (MRN with 


with / 250 yum intensity of the SPIRE 250 pm band. We used this 
higher resolution Herschel column density map for our further 
analysis of Lupus 1. 


3.1.3. Column density map from the LABOCA data 

Lupus I is sufficiently far away from ionizing UV sources to ex¬ 
clude significant amounts of free-free emission. We may there¬ 
fore assume the sub-mm fluxes to be entirely due to thermal dust 
emission. To determine whether the cloud is optically thin in 
this wavelength regime one can look at the peak intensity in the 
LABOCA map, which is = 1.37 Jy/beam. The correspond¬ 
ing optical depth can be calculated via formula (3) of [Schuller 
[etaL] ( [20091 ). 


^870yum 


-In 



ly 

Q By{T^) 


(4) 


where Q is the beam solid angle. In the case of Lupus I this yields 
values of Tgjojum ^ 0-01 foi* temperatures > 10 K. Hence 
the cloud is clearly optically thin in this wavelength regime. 
Therefore, the 870 pm intensities are directly proportional to the 
column densities of the interstellar dust and the line-of-sight ex¬ 
tinction. They can be converted to the beam-averaged hydrogen 
molecule column density via equation 

We again assume a gas-to-dust mass ratio ofR = 100 and 
for consistency extrapolate our dust model used for the Herschel 
analysis to the LABOCA wavelength, which yields a value of 
^ 870 yum = 1.32cm^g“^ Because the assumption of a constant 
dust temperature throughout the cloud is not valid, we used the 
temperature map from the Herschel analysis to calculate the 
black body function at each pixel of the map. In this way one 
obtains a column density map at the original resolution of the 
LABOCA map (i.e. EWHMgvo = 21.2"). 


3.1.4. Column density and temperature map from the Planck 
data 


^ Mathis-Rumpl- Nordsieck size distribution of interstellar grains 
( [Mathis et al.|1977| ) 


The thermal emission of interstellar dust over the whole sky 
was captured by the HEI-instrument of Planck at its 6 available 
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Fig. 2. Left: Temperature map of Lupus I from the Herschel SPIRE SED fit. Right: Histogram of the dust temperature for the whole 
cloud (purple), the northern (orange), and the center-south part (green). 


wavelengths between 350 jim and 3 mm. Together with the IRAS 
lOOyum data the emission can be well modeled by a modified 
black body. The details of the model and the fitting procedure 
are described in |Planck Collaboration et al. ( 2014 ). The resulting 
maps of temperature and optical depth (at 850 /im) can be down¬ 
loaded from the Planck archiv^ We extracted the part covering 
Lupus I from both maps to get gnomonic projected maps. The 
column density map was then computed from the optical depth 
map using the relation 


5.8 X 10^^ cm ^ ^ 

A^H 2 =--X 0.497 x 10 T 353 GHZ 


(5) 


suggested in [Planck Collaboration et al. ( 2014 ) for regions 
with higher column density than the diffuse ISM, i.e. molecu¬ 
lar clouds. 


3.2. Column density PDFs 

Probability distribution functions (PDEs hereafter) of the col¬ 
umn density are a widely used way to characterize the evolu¬ 
tion and state of molecular clouds. Both simulators (see e.g. 
Eederrath & Klessen||2012| [Ward et al.||2014[ [Girichidis et al.j 
2014|) and observers ( s ee e.g. Kainulainen et al. 2009| 20TTf 
Palmeirim et al.||2013HTremblin et al.||2014t [Schneider et al.j 
2015t Roccatagliata et al. 2015| ) use~~ttiem in their studies and 
they are a good tool to compare both simulations and observa¬ 
tions with each other. 

As Lupus I is at near distance to us and lies at a high Galactic 
latitude, far away from the Galactic plane, we expect the con¬ 
tamination of the map by overlaying foreground or background 


^ http://irsa.ipac.caltech.edu/data/Planck/release_ 
1 /al 1 - sky - map s/map s/HFI_CoinpMap_TherinalDust Mode 1_ 
2048_R1.20.fits 


emission to be small. We normalize the PDEs to the dimension¬ 
less 77 = ln(AH 2 /<A^H 2 )) which gives the opportunity to compare 
different regions of different column density, as well as the same 
cloud, but observed with dif ferent instruments. Eor s imilar obser¬ 
vations of the Orion clouds [Schneider et al.[ ( |2013| ) showed that 
the effect of varying the resolution of the maps within 18" and 
36" did not affect the main features of the PDEs (shape, width, 
etc.) significantly. The Lupus I cloud is much closer, yielding 
a much finer physical resolution and the tests we performed on 
our maps smoothing them with a grid of Gaussians with differ¬ 
ent EWHM and rebinning them to several different pixel sizes 
confirmed their findings. If the difference, however, is as big as 
e.g. the one in angular resolution between Herschel and Planck, 
it is not longer possible to maintain the features of the PDE un¬ 
changed when degrading the Herschel map by a factor of 8 
(from 500yum) to 17 (from 250yum). But in our case both 
the angular resolution and the pixel size of the Herschel and the 
LAB OCA column density map are almost the same. Therefore, 
both maps have sufficiently large numbers of pixels to have a 
high statistic. 

The basic shape of the PDE in the low column density regime 
can be described by a lognormal distribution. In more compli¬ 
cated cases the PDE might also be the combination of two log- 
normals: 


p{rf) = 



ZcTj 


I . ^2 

+ I - 

^27icr\ 


-in - Pif 

2(t\ 


( 6 ) 


where 6/ is the norm, fii and ct/ are the mean logarithmic column 
density and dispersion of each lognormal. In the high column 
density regime one often finds a deviation from the lognormal. 
This can be modeled with a power-law of slope 5- which is equiv¬ 
alent to the slope of a spherical density profile p(r) cx r~^ with 
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a 


|20T3l i. 


2/^+1 (see [Federrath & Klessen|2013[ [Schneider et al.| 


Using a least-square method (Levenberg Marquardt algo¬ 
rithm; Poisson weighting^ we derived the characteristic values 
of the distributions by fitting the PDFs with either one lognormal 
to the distribution around the single peak and where applicable 
a power-law to the high density tail, or with two lognormals to 
the distributions around the first and second peak, respectively, 
and a possible power-law tail. We checked the robustness of the 
fits by performing it on four different binsizes for the histograms 
(0.05, 0.1, 0.15, 0.2) and looking at the variation of the resulting 
fit parameters which was of the order 2% - 8%. We took the fi- 
nal parame ters from the fits for a binsize of 0.1 (see |Schneider| 
et al.||2015 ) and conservatively adopted an error of 10% for all 


fit parameters. The results are summarized in Table 


4. Results 

4.1. Column density and temperature maps 

The Herschel column density map of Lupus I is shown in 
Fig. |3(a)[ The average column density is = 1.44 x 

10^^ cm“^. The map reveals two distinct regions within the fila¬ 
ment (see also Fig.[^. The north where the column densities are 
lowest (iNH 2 6 = 1.10x10^^ cm“^) and just one dense core can be 
seen ( |Rygl et ^|2013| ) and the center-south part where several 
cores with high column densities above 10^^ cm“^ are found and 
the average column density = 1.72 x 10^^ cm“^) is about 

60% higher than in the northern part. 

Comparing our map to the one derived by |Rygl et al. ( 2013| ) 
we note that our values are lower by a factor of about 2-2.5. This 
can be explained by the use of a different dust model (opacity 
and dust spectral index) and the inclusion of the PACS 160yum 
band in their work. The extinction maps created from 2MASS 
and Spitzer data (see [Chapman et al. [2007 [ [Lombardi et al.|2008[ 
[Merm et al.|2008| ) have a resolution that is lower by a factor of 
~ 7-30. Nevertheless, between those maps and our column 
density map is a clear resemblance. 

The column d ensity map of Lupus I we obtained from 
LABOCA (see Fig. [3(b)| ) is dominated by the dozen cores in the 
center-south part and the denser dust in the central part of the fil¬ 
ament. The average column density is = 1.47 x 10^^ cm“^. 

The Planck column density map of Lupus I with a resolution 
of 5' reveals the basic structure of the cloud which agrees well 
with the column density maps from LABOCA and Herschel. 
This is consistent with what is shown by the other two maps 
when smoothed to the resolution of Planck and thus smoothing 
out the highest column density peaks. 

The temperature map obtained from the Herschel SED fit 
(see Fig. [2(a)| ) shows an anti-correlation of the temperature with 
the column density. The densest parts are the coldest and the less 
dense the material the warmer it becomes. Also one sees again 
a difference between the northern and the center-south part of 
the cloud. In the north dust temperatures between 15 - 23 K 
with a mean and median of 19 K are found. But only in the 
dense pre-stellar core and slightly north of it the temperature 
drops down to 15 K. The maximum temperature in the center- 
south part is 22 K at the edges of the filament. The inner part of 
the filament has an average temperature of 17 K. But it drops 
down to even 11 K in the densest cores. From the histogram 
of the dust temperatures (see Fig. 2(b)[ ) one can see that most 
of the dust in the north (orange histogram) ranges between 
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18 - 20 K whereas in the center-south part (green) it is between 
16 - 18 K. This means that the dust in the center-south is on 
average 2 K colder than in the north. 



Fig. 4. 3D-surface plot of the Herschel 250yum column density 
map. One can clearly identify the cores in the center-south of 
Lupus I and see the on average lower column density in the 
northern part of the cloud. The red dashed boxes mark the two 
regions we distinguish in this work. 


4.2. Column density PDFs 

We derived PDFs of the Lupus I cloud from the Herschel and 
LABOCA column density maps for the two distinct parts of the 
cloud (north and center-south; see Fig. |^, as well as for the 
whole cloud. They are shown in Fig. For a correct interpre¬ 
tation of the PDFs one has to consider the completeness limit 
of the underlying column density map. Here we adopted the 
lowest closed contour as such a limit which is 1.1 x 10^^ and 
2x10^^ cm“^ for the Herschel and LABOCA map, respectively. 
It is marked by the vertical dashed line in the plots of Fig.[^ 
Looking at the entire cloud the Herschel PDF is very com¬ 
plex with clear deviations from a simple lognormal distribu¬ 
tion which would be expected for a cloud that is dominated by 
isothermal, hydrodynamic turbulence (e.g. [Klessen et al.|2000] ). 
The distribution shows two peaks in the low column density 
regime and a power-law tail in the high density end. The first 
peak at NH 2 ,peaki = 8.66 x 10^^ cm"^ falls below our complete¬ 
ness limit, so it might be just refiecting the drop of observa¬ 
tional sensitivity and will not be used for interpretation. But the 
distribution above the limit can be well represented by a fit of 
two lognormals around the first and second peak which is at 
NH2,peak2 = 2.57 X 10^^ cm“^. The width of the first compo¬ 
nent is (Ti = 0.56, but could also be broader due to the possi- 
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Table 2. Results of the fits to the Herschel PDFs of Fig. 


Region 

o-i 

Nh 2 ,peakl 

[xlO^^ cm"^] 

0-2 

NH2,peak2 

[xl02‘cm-2] 

< Nh 2 > 
[xlO^^ cm"^] 

NH2,dev 

[xl02‘cm-2] 

s 

a 

Whole cloud 

0.56 

0.866 

0.26 

2.57 

1.44 

3.60 

-2.55 

1.78 

Center-south 

0.46 

0.934 

0.35 

2.71 

1.72 

4.50 

-2.43 

1.82 

north 

0.55 

0.794 

0.18 

2.23 

1.10 

- 

- 

- 



Fig. 3. Column density maps of Lupus I from (a) the Herschel 250/im map, (b) LABOCA, and (c) Planck. 


ble underestimate of the column densities below the complete¬ 
ness limit. Nevertheless, it is more than twice as broad as the 
width of the second lognormal ((T 2 = 0.26) and larger than in 
other nearby clo uds like Maddalena, A uriga ( [Schneider et al.| 
|2Q13| ) or Aquila ( [Schneider et al.|[2QT5] ). This can be a sign of 
broadening by turbulence and external compressive forcing (see 


Federrath et al.|201^[Tremblin et al.|20T^[Federrath & Klessen[ 


2013[ ). The power-law tail follows a slope of s = -2.55 or cor¬ 
responding a = 1.78 for a spherical density distribution. The 
deviation from the lognormal into the power-law occurs near 
NH 2 ,dev = 3.6 X 10^^ cm“^ which is slightly below the recent 


values by [Schneider et al.[ ( 2015| ) who found the transition into 
the power-law tail to be at Ay ~ 4 - 5 ma^ in their investiga¬ 
tion of four low-mass and high-mass star forming regions. But it 
is higher than the result of [Roccatagliata et~SI| ( [2Q15[ ) who find 
the deviation to occur at Ay ~ 1.2-2 mag in the Serpens Core 
region. 

From the Herschel data the center-south part of the cloud 
shows a PDF very similar to that of the entire cloud. Here again 
two peaks are present and the distribution above the complete¬ 
ness limit can be fitted by a double lognormal and a power-law 
tail. The width of the first fitted lognormal is narrower, the width 


N(H 2 )/Av = 0.94 X 10^^ cm ^mag ^ ( Bohlin et al.|l978 1. 


of the second broader than those of the entire cloud with val¬ 
ues of (Ti = 0.46 and (T 2 = 0.35. The slope of the power-law 
is slightly flatter with s = -2.43, but leading to a comparable 
exponent a = 1.82. 

Also in the northern part of the cloud the double peak log¬ 
normal profile of the PDF from Herschel is found again. But 
now there is no excess above the second lognormal at higher 
densities. The peak position of the first lognormal is lower 
(NH 2 ,peaki = 7.94 X 10^^ cm“^) compared to the center-south and 
the whole cloud, but with similar width (cri = 0.55). The peak 
position of the second lognormal (NH 2 ,peak 2 = 2.23 x 10^^ cm“^) 
is lower than that of the whole cloud and its width ( 0-2 = 0.18) 
is approximately half as broad. 

Although it is not possible to strictly constrain the shape of 
the PDF below the completeness limit, one can say that the dis¬ 
tribution above shows not just two components (lognormal and 
a power-law) but at least three because of the second peak that 
appears before the deviation into the power-law. Such a shape 
could be attributed to an initial turbulent cloud which later was 
compressed by an external driving agent. Such an agent can ei¬ 
ther be an ionization front, an expanding shell driven by winds 
or/and supernovae or colliding flows. The initial lognormal form 
of the PDF of the cloud develops a second component caused by 
the dense compressed gas. This behavior was found in simula- 
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tions by Tremblin et al.| ( [20T2 ) of an initial turbulent medium that 
is ionized and heated by an ionization source. They concluded 
that a double peak is present when the ionized gas pressure dom¬ 


inates over t he turbulent ram pressure of the cloud. [Matsumoto 
|et al.| ( |MT5] ) studied the evolution of turbulent molecular clouds 
swept by a colliding flow. They found that the PDF exhibits two 
peaks in the case when the Mach numbers of the initial tur¬ 
bulence and the colliding flow are of the same order M = 10 
(model HTIOFIO; see their Fig. 11b) and the line of sight is per¬ 
pendicular to the colliding flow. Then the low column density 
peak represents the colliding flow and the higher column den¬ 
sity peak the sheet cloud, respectively. Qbservationally, double 


peak PDFs have been reported an d stu died recently by |Schneider 
|et al. I ( [20T^ , [Harvey et al.| ( |2Q13| ), and |Tremblin et al.| ( |2Q14| ) for 
several nearby clouds exposed to an ionization source. Areas of 
the cloud close to the ionization front indeed showed the pre¬ 
dicted second peak in the PDF due to the compression induced 
by the expansion of the ionized gas into the molecular cloud. In 
the case of Lupus I the source of the compression is very likely 
the expanding HI shell around the USco sub-group of Sco-Cen 
(possibly together with a supernova that exploded within USco) 
and the wind bubble of the remaining B-stars of UCL pushing 
from the eastern and western side of the cloud, respectively. This 
scenario will be discussed in more detail in Sec.[5l 

Lognormal PDFs indic ate shock waves ([Kevlahan & Pudritz 


20091 or turbulence (e.g. 
et al. 2010| ). [Tremblin et a. 


Vazquez-Semadeni 1994[ Federrath 
.|( |2012| ) have shown in a simulation 
of Stroemgren spheres advancing into turbulent regions that the 
region compressed by the ionization front also has a lognormal 
density PDF, but shifted to higher density due to the compres¬ 
sion. A PDF of a region that includes compressed and undis¬ 
turbed parts of the cloud will thus show a double peak PDF. 

The power-law tail that is seen in our PDFs could be ex¬ 
plained by active star formation and the transition to a gravity 
dominated density regime represented by the star forming cores. 
This can be shown by comparing the PDFs of the two regions in 
the north and the center-south. In the north where there is almost 
no star formation and just one pre-stellar core and 20 unbound 
cores can be found ( jRygl et al.|2013| ), the PDF shows no power- 
law tail in the high density regime. In the center-south instead 
where almost all the star formation activity takes place with 
plenty of very young bound cores, the PDF shows the power- 
law tail very clearly. Also numerical studies have shown that the 
PDF for an actively star forming region develops a clear devia¬ 
tion from the lognormal in form of a power-law tail ( jBallesteros- 
jParedes et al.|201ll|Girichidis et al.|2()T4] ). 

With LAB OCA the PDF of Lupus I can be modeled as one 
lognormal with a power-law tail at high densities (NH 2 ,dev = 
6.0 X 10^^ cm“^). The peak of the lognormal distribution lies at 
NH2,peak = 1.43 X 10^^ cm"^ which is below the completeness 
limit of the LAB OCA column density map. The width of the 
lognormal is broad (cr = 0.56), but uncertain due to the peak 
position and thus possibly even broader. The power-law slope 
(s = -2.66, a = 1.75) is very similar to the one of the Herschel 
PDF indicating the star formation activity in the cloud. 

The PDF of the center-south region shows the same be¬ 
havior as the one for the whole cloud. The lognormal part is 
slightly broader (cr = 0.60), but the power-law slope {s = -2.59, 
a = 1.77) and the deviation point (6.0 x 10^^ cm“^) are almost 
the same compared to the whole cloud. In the northern part only 
the now narrower lognormal is seen instead (cr = 0.52), indicat¬ 
ing the lack of star formation activity and the dominance by the 
turbulence induced trough the compression. In all three cases the 
positions of the peak are very close to each other. 


As was already seen in the column density map from 
LAB OCA, its PDF shows only the dense cores and the intercore 
medium. The instrument is not sensitive to all material that can 
still be probed with Herschel. With LAB OCA we thus sample 
more the second lognormal of the Herschel PDF and the power- 
law tail. 


4.3. Core distribution 


For the core analysis in the LAB OCA map we used the 
Clump find package ( [Williams et al.||1994| ). It decomposes the 
emission of the map into a set of clumps or cores contouring 
the data at given threshold levels. The results can be found in 
Tab. and the distribution of the cores (represented by the white 
crosses) is shown in Fig.[2 The algorithm identifled 15 different 
cores with masses between 0.07 and 1.71 M©. Their total mass 
is 8.37 Mq and their total flux is 25.39 Jy. This corresponds to 
^ 5% of the total mass of Lupus I (see Sec. [4.4[ ) and 5% of 
the total flux of the LAB OCA map, respectively. For the com¬ 
putation of the core masses, we derived their temperatures from 
the Herschel SPIRE SED fit temperature map as a mean tem¬ 
perature within the ellipse representing the core. With this tem¬ 
perature the mass of the cores was calculated following [Schuller 
[etaL] 120091 ) 


d^FyR 

v(T^)k 


where d is the distance to Lupus I, Fy the total flux of the 
core, R = 100 the dust-to-gas ratio, and Ky = 1.32cm^g“^ (see 
Sec. 3.1.3[ ). Most interestingly, all but one core are found in the 
center-south part of Lupus I. This confirms our interpretation of 
the column density maps that is the part where active star forma¬ 
tion is going. _ 

T he dist ribution of the cores from the Herschel data ( |Rygl 


et al.[[2013] ) agrees with this picture and the core distribution 


from our LAB OCA map. The northern part of Lupus I is mainly 
populated by unbound cores and just one pre-stellar core wheras 
the center-south is dominated by pre-stellar cores. Many of those 
coincide with our LAB OCA cores. However, a direct assignment 
is not easily possible since the Herschel coordinat es can only 
be estimated from the source images in the paper of [Rygl et al. 


( [2013[ ) and in some cases several Herschel sources seem to be 
on the position of one LAB OCA source or vice-versa. 

Benedettini et al.[ found eight dense cores in Lupus I 


using high-density molecular tracers at 3 and 12 mm with the 
Mopra telescope (red diamonds in Figurej^. Seven of those have 
one or more counterparts in our LAB OCA map. The matches are 
given in Table They classify five of their cores as very young 
protostars or pre-stellar cores (Lupl C1-C3, C5, and C8) and 
the remaining three also as very likely to be protostellar or pre- 
stellar. For three of their cores (Lupl C4, C6, and Cl) they cal¬ 
culated a kinetic temperature of 12 K. This agrees within 20% 
with the core temperatures derived from the Herschel SPIRE 
SED fit temperature map. 

From the 17 young stellar ob jects found in Lupus I by the 
Spitzer c2d near-infrared survey ( [Chapman et al.[[2007} Merm 


et al.[[2008 ) 11 lie within the boundaries of the LAB OCA map 


(blue boxes in Figure [^. But only one object (IRAS 15398- 
3359) clearly matches one of our LABOCA cores (#1). This is 
the long known low-mass class 0 protostar IRAS 15398-3359 
(inside the B228 core) which has a molecular outflow (see re¬ 
cently e.g. [Oya et al.|2014[ Dunham et al.|2014] ). Two other ob¬ 
jects of the c2d survey are close, but offset by 0.5' to the center 
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of our cores #3 and #5, respectively. Nevertheless, each of those 
two objects still lies within the boundaries of the corresponding 
ellipse representing the LAB OCA core (see Table |^. Most of 
the LAB OCA cores are potentially at a very early evolutionary 
stage, i.e. without a protostar inside to heat it and eventually de¬ 
stroy the surrounding dust envelope to allow near-infrared radi¬ 
ation to escape. In fact the only cores with a 70/im counterpart, 
which is a good proxy for a protostar inside a core, are num¬ 
ber #1 and #3. Therefore, we do not expect to find many Spitzer 
counterparts. 



Fig. 6. Herschel column density map of the center-south region 
with contour levels of 1.1, 4, 7, and 10 x 10^^ cm“^. The differ¬ 
ent symbols mark the center positions of the cores found with 
LAB OCA and of objects found in other studies as discussed in 
Section l43l 


Vilas-Boas et al. ( 20Q0| found 15(14) condensations in C^^O 


(^^CO) with the 15 m SEST telescope. Eight cores lie within our 
LABOCA map (orange crosses in Eigurej^ and of those three 
(Lu7, LulO, B228) coincide with LABOCA detected cores (#12, 
#5, and #1). Their C^^O excitation temperatures are 9, 11, and 
10 K, respectively. This is 20 - 40% lower than the tempera¬ 
tures from the Herschel map. Eor cores Lu7 and LulO they cal¬ 
culated a mass from C^^O of 11.7 and 3.1 M©, respectively. This 
is a factor of 20 and 7 higher than the LABOCA masses. 
But the sizes of their condensations were on average larger than 
at least three times their beam size of 48". 


Eigure shows the Herschel column density map of the 
center-south region with contour levels of 1.1, 4, 7, and 10 x 


10^^ cm“^. Overploted are the LABOCA cores found in this 
study (green cir cles), as well as the abov e mentioned cor es and 
YSOs foun d by|Benedettini et al. ( |2012| ) (red diamonds), Vilas- 
Boas et aL] ( |2000| ) (orange crosses), and the Spitzer c2d survey 
(blue boxes), respectively. 

Our findings confirm that Lupus I harbors a population of 
very young cores that are forming stars just now and probably 
have ages below 1 Myr. This could support the idea of an ex¬ 
ternal shock agent, like the USco shell, sweaping up the cloud 
and triggering the simultaneous formation of new stars within 
the cloud. 


4.4. Total mass estimates of Lupus I 

To derive the total mass of the cloud from the three different 
column density maps we defined a polygon around the filament 
(delineated in Eig.[2on the LABOCA map) to derive the total 
mass always in the same area 1 deg^). The total gas and dust 
mass was then calculated via the formula 


Lupus 


I - ^ Nh2 PII 2 Ap 


( 8 ) 


with Ap the area of a pixel in cm^. As the common lower level we 
chose Nh 2 >10^^ cm“^ which corresponds to Ay > 1 mag. The 
resulting total mass for Lupus I is j 171 M©, 174 M©, 

and 164M© for the Planck, Herschel, and LABOCA data. 
This means that the total masses calculated from the three data 
sets agree with each other. Comparing the total mass of the cloud 
to the total mass in cores from LABOCA 8 M©) one sees that 
only about 5% of the mass is concentrated in the densest conden¬ 
sations. Most of the dust and gas is in more diffuse components. 

The most recent literatur e value for the tota l mass of Lupus I 
is from the Herschel data by Rygl et al. ( 2013| ). They calculated 
the total mass in a much bigger area than this work (^^ 4.5 deg^ 
compared to our 1 deg^) and used a different method to cre¬ 
ate their column density map as mentioned already in Sec. |4.1[ 
Therefore, a difference of about a factor 4-5 arises between 
their value of M = 830 M© (for Ay > 2 mag) and our finding. 

Other literature values cover a wide range of total masses 
for Lupus I, depending on the tracer used and the size of the 
area that was considered. Erom the Spitzer c2d near-infrared ex¬ 
tinction maps |Merm et~n^ (|2008|) determined a total mass of 


479 M© for Ay > 3. Various CO measurements (e .g. |Tachihara 
et al.| 1996[ Kara et aL]|1999t [Tothill et~aL]|2009| ) yielded val¬ 
ues of 280 - 880 M©. Direct comparisons with our values are 
not always possible since all maps cover different parts of the 
Lupus I cloud complex and the material is not homogeneously 
distributed to allow scaling the mass with the area. But we note 
that our values agree with most of the literature values within a 
factor of 2-3 which is expected considering the uncertainties in 
the choice of the dust model, the dust-to-gas ratio and the CO- 
to-H 2 conversion factors. 


5. The surroundings of Lupus I and the interaction 
with USco and UCL 

Best suited for looking at the dust surroundings of Lupus I are 
the Planck data (see Eig. [7] and Eig. [^. These observations at 
350, 550, and 850yum cover the whole sky at a resolution of 
5'. Lupus I lies on the eastern edge (Galactic coordinates) of 
a ring-like dust ridge (labeled in Eig. that extends from about 
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Table 3. Parameters of the 15 cores detected in the LABOCA map with Clumpfind (see Fig.[^. The beam size (HPBW) and area 
used to convert intensities into fluxes is 21.2" and 509 arcsec^, respectively. The total mass of all cores is 8 M©. The 4th and 5th 
column give the FWHM of the ellipses representing the cores. The 6th, 7th, and 8th column give the core’s peak intensity, total flux, 
and mass, respectively. The core temperature shown in the 9th column was derived from the Herschel SPIRE SED fit temperature 
map as a mean temperature within the ellipse representing the core. The last column gives the matches to dense cores found by 
[Benedettini et al.| ( |2012| ). 


# 

Ra 

Dec 

FWHM^ 

["] 

FWHMy 

["] 

Imax 

[Jy/beam] 

Ftot 

[Jy] 

Mcore 

[Me] 

Tcore 

[K] 

Counterpart 

1 

15:43:01.68 

-34:09:08.9 

58.22 

44.38 

1.37 

5.61 

1.71 

14.0 

Lupl C4 

2 

15:44:59.85 

-34:17:08.8 

62.51 

44.09 

0.51 

3.97 

1.59 

12.0 

Lupl C6 

3 

15:45:13.58 

-34:17:08.8 

46.56 

33.45 

0.33 

1.48 

0.40 

15.0 

Lupl Cl 

4 

15:45:16.02 

-34:16:56.6 

35.31 

34.95 

0.31 

1.00 

0.27 

15.0 

Lupl C7 

5 

15:45:03.80 

-34:17:57.3 

57.22 

42.24 

0.28 

1.41 

0.46 

13.5 

Lupl C6 

6 

15:42:52.86 

-34:08:02.1 

96.42 

49.94 

0.28 

3.85 

1.34 

13.0 

Lupl C3 

7 

15:42:43.57 

-34:08:14.2 

70.62 

44.38 

0.25 

1.64 

0.57 

13.0 

Lupl C3 

8 

15:42:44.55 

-34:08:32.5 

30.55 

20.85 

0.25 

0.48 

0.17 

13.0 

Lupl C3 

9 

15:42:44.06 

-34:08:44.6 

36.82 

19.26 

0.24 

0.49 

0.17 

13.0 

Lupl C3 

10 

15:42:50.42 

-34:08:38.5 

53.28 

30.34 

0.22 

1.05 

0.36 

13.0 

Lupl C3 

11 

15:45:25.10 

-34:24:01.8 

42.46 

69.57 

0.22 

1.63 

0.53 

13.5 

Lupl C8 

12 

15:42:23.02 

-34:09:33.2 

68.97 

46.08 

0.21 

1.72 

0.53 

14.0 

Lupl Cl 

13 

15:39:09.92 

-33:25:30.9 

27.24 

27.70 

0.20 

0.34 

0.07 

17.0 

- 

14 

15:45:55.75 

-34:29:23.7 

28.10 

26.28 

0.20 

0.29 

0.08 

14.5 

- 

15 

15:42:47.90 

-33:53:15.3 

43.80 

19.59 

0.19 

0.43 

0.12 

15.0 

Lupl C2 


/? = +10° to /? = +25° in Galactic latitude with a center at about 
/ = +345°, b = 17.5°. The ridge is ~ 5° wide which corresponds 
to ~ 13 pc at the distance of Lupus 1. Besides Lupus I it consists 
of several small molecular clouds extending north of Lupus I and 
then bending towards the east connecting with the p Ophiuchus 
molecular cloud on the opposite site of Lupus I. Further west of 
Lupus I no dust emission is seen in the Planck maps. The same 
is true for the inside of the dust ridge. Between p Ophiuchus 
and Lupus I a roundish dust void is seen. However, these two 
dust voids on either side of Lupus I are fllled with hot X-ray gas 
which can be seen with ROSAT (left in Fig. [7]). Fig.j^shows the 
dust emission seen by Planck in 350//m. Overlayed in white are 
the contours of the ROSAT diffuse X-ray emission at 3/4 keV[^ 
Inside USco the contours follow the edge of the dust ridge indi¬ 
cating that the hot X-ray gas might be in contact with the cold 
dust. On the western side of Lupus I the contours mark the out¬ 
line of the second roundish X-ray emission seen in ROSAT. The 
cyan dots mark the remaining early B-type stars of UCL. Both 
contours seem to wrap around Lupus I what might be a sign that 
the cloud is embedded in hot ISM. 

These dust voids and the observed X-ray gas might be ex¬ 
plained by the cumulative feedback of the massive stars in the 
USco and UCL sub-groups of Sco-Cen. Their creation has been 
interpreted in a scenario of propagating molecular cloud for¬ 
mation and triggered cloud collapse and star formation taking 
place within the last 17 Myr in the UCL and USco sub-groups of 
Scorp ius-Centaurus (see |de Geus||1992[ Preibisch & Zinnecker| 
|2007| ). The expansion of the UCL HI shell that started ~ 10 Myr 
ago and was driven by winds of the massive stars and super¬ 
novae explosions has probably cleared out almost all the dust 
and molecular material west of Lupus 1. But the observed X- 
ray gas on that side is probably not related to the supernova 
explosions, because the X-ray luminosity of supernova-heated 
superbubbles dims on a timescale of typically less than a mil- 


lion years via expan sion losses and mixing with entrained gas 
( Krause et al.||2014| ). The shell today has a radius of ~ 110 pc 
(as seen from HI data; |de Geus||1992|), general ly consistent 
(comp are, e.g. [Baumgartner &^eitschwerdt|20 13} [Krause et al.| 


2013) with the inferred age of the stellar group, about 17 million 


^ The very strong X-ray emission seen near the position of the B2IV 
star south-west of Lupus I (/cCen; / = 326.872°, h - 14.754°) is neither 
related to the star nor to Sco-Cen, but most likely caused by the Quasar 
[VV2006] J150255.2-415430. 


years. Therefore, this X-ray gas is probably currently heated by 
the remaining B-stars of UCL. Six of those stars south-west of 
Lupus I (325° < I < 335° and 10° < b < 15°; see Fig.[|( lie at 
positions that favor them as likely being the sources of this wind 
bubble. 

The USco HI shell has started its expansion ~ 5 Myr ago 
powered by the winds of the OB-stars and quite possibly a recent 
supernova explosion, about 1.5 Myr ago, as suggested by the de¬ 
tection of the 1.8 MeV gamma ray line towards the USco sub- 
group and the detection of the pulsar PSR J1932+1059 ( [Diehl 
jet al.|2010| and references therein). What now forms the before 
mentioned dust ridge is probably the remaining material of the 
USco parental molecular cloud that was swept up by the advanc¬ 
ing USco shell leaving the dust void fllled with X-ray gas. From 
simulations, we expect large scale oscillations of the hot. X-ray 
emitting gas ( [Krause et al.||2014| ), because the energy source is 
never exactly symmetric. 

It seems like the UCL wind bubble could be colliding with 
the USco HI shell right at the position of Lupus I squeezing it 
in-between. This wind bubble might have provided a counter¬ 
pressure to the expanding USco shell and thus favored this po¬ 
sition for an additional compression of the shell material. In this 
way a new molecular cloud could have been created there and it 
might explain why we do not see more very young star forming 
clouds (except p Oph) distributed within the wall of the USco 
shell. With its distance of 150 pc the cloud also is neither in 
the foreground of t he B-stars of UCL that have an average dis¬ 
tance of ~ 140 pc (de Zeeuw et al.[[1999[ ) nor in front of USco 
(~ 145 pc; Preibisch et al.[2002| ). To rule out shadow effects 
completely, however, a more comprehensive analysis of the X- 
ray data would be required. 

From our preliminary analysis of the GAS S HI data 
( [McClure-Griffiths et al.|2009l[Kalberla et al.|2010[ ) we see that 
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Lupus I is indeed embedded within the expanding shell around 
US CO (Kroell et aL, in prep.). The shell has an expansion veloc¬ 
ity of 7kms“^ and a thickness of 6 pc. We estimate the 
current outer radius to be 36 pc (both inner and outer radius of 
the USco shell from our HI model are delineated in Fig.|^. 

These findings complement previous suggestions of an in¬ 
teraction of Lupus I (and p Oph) with the two sub-gro ups of 
Sco-Cen that were originally proposed by |de Geus ( 1992|). From 
the spatial and velocity structure of observations |Tachihara 


let al.|(200T] ) also found evidence for an interaction. |Tothill et al. 
( 2009 ) in their analysis of and CO(4-3) found enhanced 
line widths at the western end of Lupus I and a velocity gradient 
across the filament, i.e. in the direction of the USco shell expan¬ 
sion, of ~ 1 kms"^ which they conclude to be consistent with a 
dynamical interaction between Lupus I and the USco HI shell. 
Matthews et al.| ( [MT4] ) have found that the large-scale magnetic 


field is perpendicular to the Lupus I filament, i.e. points in the 
direction of the USco shell expansion. This might have favored 
the accumulation of cold, dense atomic gas along the field lines 
and promote faster molecule fo rmation ( jHartmann et al.|[^01[ 
Vazquez-Semadeni et al.|20lT ). 


All these are hints towards the scenario that Lupus I was af¬ 
fected by large-scale external compressing forces coming from 
the expansion of the USco HI shell and the UCL wind bubble. 
This might explain its position, orientation and elongated shape, 
the appearance of a double peak in the PDF and the large pop¬ 
ulation of very young pre-stellar cores that are seen with both 
LAB OCA and Herschel. 


6. Summary and Conclusions 

From our LAB OCA observations of Lupus I as well as from 
archival Herschel and Planck data we created column density 
maps of the cloud. In addition we calculated a temperature map 
from an SED fit using the three SPIRE bands of Herschel. All 
maps suggest that the cloud can be divided into two distinct re¬ 
gions. The northern part that has on average lower densities and 
higher temperatures as well as no active star formation and the 
center-south part with dozens of pre-stellar cores where density 
and temperature reach their maximum and minimum, respec¬ 
tively. 

From the LAB OCA and Herschel maps we derived column 
density PDFs for the entire cloud as well as for the two above 
mentioned regions separately. The Herschel PDF of Lupus I 
showed a double peak profile with a power-law tail. The power- 
law can be attributed to the star formation activity in the center- 
south part of the cloud since it disappears in the PDF of the 
northern part. However, the double peak profile is conserved 
throughout the cloud and possibly arises due to the large-scale 
compressions from the cumulative massive star feedback of the 
Sco-Cen sub-groups. Such a sign for compression in the PDF 
was previously found in both observations and simulations of ad¬ 
vancing ionization fronts, also supported by simulations of col¬ 
liding turbulent flows. With LAB OCA we probe the denser parts 
of Lupus I and find one lognormal and the power-law tail for the 
whole cloud and the center-south. The PDF of the northern part 
shows only the lognormal behavior. 

The distribution of the 15 cores that were found in our 
LAB OCA map confirms that only the center-south part of 
Lupus I is actively forming stars whereas the north is quiescent. 
The cores have masses between 0.07 and 1.71 M© and a total 
mass of 8 M©. As the cloud has a total gas and dust mass of 
^164 Mq (from LAB OCA for Nh 2 >10^^ cm“^) this means that 
~ 5% of the mass is in cores. All of those cores are pre-stellar 


or protostellar confir ming that we are witnessing a large star for - 
mation event (see e.g. |Benedettini et al.|2012t|Rygl et al.[2013| ). 
The total mass from Herschel and Planck for Nh 2 > 10"^^ cm“^ 
was 174 M© and 171 M©, respectively. 

We argue that the main driving agents in the formation pro¬ 
cess of Lupus I are the advancing USco HI shell, which most 
likely contains the cloud, and the UCL wind bubble interacting 
with the USco shell and hence squeezing Lupus I in-between 
these two shells in a complex manner. The age of the pop¬ 
ulation of young stellar objects suggests a compression event 
~ 1 - 2 Myr ago, which could be the interaction between the 
two shells and possib ly a supernova exp losion as suggested by 
the gamma ray data ( |Diehl et ^|2010| ). This large-scale com¬ 
pressions might be the reason for the position, orientation and 
elongated shape of Lupus I, as well as the double peak PDFs 
and the population of very young pre-stellar cores we found. 

In future work we will analyze our newly performed high- 
resolution APEX CO line observations to study the kinematics of 
the Lupus I cloud. We will also compare our observational find¬ 
ings with dedicated numerical simulations and search for signs 
in the theoretical PDF that will allow us to distinguish between 
different scenarios of molecular cloud formation. 
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Fig. 5. Column density PDFs of Lupus I and their model fits for the north and the center-south part of the cloud, as well as for 
the entire cloud. The error bars show the ^^N uncertainties. In case of Herschel the two yellow and purple dashed lines show the 
fits of the two lognormals to the distribution around the first and second peak, respectively. Their superposition is represented by 
the solid green line. For LAB OCA the solid green line shows the lognormal fit to the distribution around the peak. For all cases 
the straight solid red line shows the power-law fit to the high-density tail of the PDF, if applicable. In each plot we give the value 
of the dispersion(s) of the fitted lognormal(s) ct/, the position(s) of the peak(s) of the lognormal(s) NH 2 ,peaki 7 the slope of the fitted 
power-law tail s, the corresponding slope of an equivalent spherical density distribution a, the mean column density and the 

deviation point from the lognormal to the power-law tail NH 2 ,dev The vertical dashed line marks the incompleteness limit (lowest 
closed contour) which is 1.1 x 10^^ and 2 x 10^^ cm“^ for the Herschel and LABOCA map, respectively. 
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Fig. 7. The surrounding of Lupus I in X-rays (ROSAT 3/4 keV, left) and dust emission {Planck 850yum, right). There are two X-ray 
bubbles on either side of the Lupus I ridge and both Lupus I and p Ophiuchus as well as the little cloudlets north of them seen in 
dust emission have X-ray shadows. The yellow dashed circles in the right image give the position of the inner and outer edge of the 
USco HI shell determined from our preliminary analysis of the GASS HI data. The center lies at I = 347°, h = 25° and the inner 
and outer radius is 12° and 15°, respectively. 



Fig. 8. The surroundings of Lupus I seen in dust emission {Planck 350/im). The white contours mark the X-ray emission (ROSAT 
0.8 keV) from 0.27 -lx 10“^ counts/s. The cyan dots mark the positions of the B1-B3 stars of UCL from the Hipparcos catalog. 
Labels are given in green. 
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